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-  Pedersen  conductivity  gradient  scale  lengths  L  =  4,  6, 10  km  perpendicular  to  the  mag¬ 
netic  field.  Independent  of  the  initial  conditions  used,  all  of  the  models  studied  become 
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with  ny  —  2—3  for  27r/ky  between  1  and  10  km.  These  quasi-final  power  law  power 
spectrum  states  are  achieved  faster  for  the  random  initial  perturbation  than  for  the 
monochromatic  initial  perturbation,  and  similarly  for  the  large  amplitude  case  within 
each  set. 
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EFFECT  OF  DIFFERENT  INITIAL  CONDITIONS  ON  THE 
EVOLUTION  OF  THE  E  X  B  GRADIENT  DRIFT  INSTABILITY  IN 
IONOSPHERIC  PLASMA  CLOUDS 

1.  INTRODUCTION 

Experimental  studies  of  artificially  injected  plasma  clouds  into  the 
ionosphere  [ Rosenberg,  1971:  Davis  et  al . ,  1974;  Baker  and  Ulwick,  1978] 
have  provided  much  information  concerning  not  only  ambient  ionospheric  con¬ 
ditions,  e .  g.,  electric  and  magnetic  fields,  but  also  the  structure  and 
morphology  of  evolving  plasma  clouds  themselves  by  means  of,  for  example, 
scintillation  and  power  spectrum  studies*  The  characteristic  initial 
steepening,  elongation,  and  striation  of  drifting  plasma  clouds  have  been 
explained  by  applying  the  linear  theory  of  the  ExB  gradient  drift  Insta¬ 
bility,  originally  developed  for  laboratory  plasmas  [simon,  1963],  to  plasma 
cloud  geometries  [Haerendel  et  al * ,  1967;  Linson  and  Workman,  1970; 

Volk  and  Haerendel,  1971;  Perkins  et  al . ,  1973].  More  recently , Chat urvedi 
and  Ossakow  (1979)  have  studied  the  nonlinear  stabilization  of  the  long 
wavelength  ExB^  gradient  drift  instability  in  ionospheric  plasma  clouds. 

Numerical  simulation  studies  f Zabusky  et  al . ,  1973;  Lloyd  and 
Haerendel ,  1973;  Goldman  et  al . ,  1974;  Doles  et  al . ,  1976:  Ossakow 

et  al . ,  1975;  197?]  of  barium  plasma  clouds  with  a  background  ionosphere 
have  reproduced  not  only  many  of  the  gross  observat ional  features  of 
plasma  cloud  evolution,  but  also  their  spatial  power  spectra  [ Scannapieco 
et  al. ,  1976]  ,  minimum  scale  size  [McDonald  et  al . ,  1978;  1980]  and 
outer  scale  size  or  correlation  length  [ Keskinen  et  al . ,  1980a]  .  In 
addition  numerical  simulations  of  the  local  ExB  gradient  drift  insta¬ 
bility  in  ionospheric  plasma  clouds  [Keskinen  et  al.,  1980b  ]  have  yielded 
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Spatial  power  spectra  and  saturated  wave  amplitudes  that  are  consistent 
with  experiment. 

However,  to  date,  the  effects  of  different  initial  conditions,  i.e., 
changing  the  initial  seed  perturbations,  on  the  spatial  and  temporal  evolu¬ 
tion  of  ionospheric  plasma  clouds  have  not  been  studied  in  detail.  The 
goal  of  the  present  work  is  to  study  the  effects  of  varying  the  initial 
conditions  on  the  evolution  of  larger  region  initially  slablike  ionospheric 
plasma  clouds  both  in  real  (x,y)  space  and  in  Fourier  (kx,  kv)  space.  In 
rection  2  we  present  the  model  equations  appropriate  for  describing  large 
plasma  clouds  in  the  F  region  ionosphere.  The  results  of  the  numerical 
solution  of  these  equations  under  several  different  sets  of  initial  con¬ 
ditions  are  given  in  sections  3  and  4.  A  summary  and  discussion  of  the 
principal  results  of  this  study  are  presented  in  Section  5. 

2.  MODEL  EQUATIONS 

For  wavelengths  much  greater  than  the  ion  gvroradius  (approximately 
10  meters  for  Ba*  ions  in  the  twilight  F  region),  the  dynamics  of  the 
plasma  cloud  and  background  ionosphere  can  be  studied  in  the  fluid  approxi¬ 
mation  (Vo* lk  and  Haerendel ,  1971;  Perkins  et  al . ,  1973,  Zabusky  et  al . , 
1973;  Ossakow  et  al . ,  1 9 / 5 ]  .  For  large  clouds  (large  magnetic  field  line 
integrated  Pedersen  conductivity  compared  with  that  of  the  background  iono¬ 
sphere),  the  cloud  interaction  with  the  background  ionosphere  (second  level) 
can  be  neglected  [  Haerendel  et  al . ,  1967  ]  ,  Furthermore,  due  to  the 
very  high  conductivity  along  the  magnetic  field  lines  (typically  ap^ a  jj 
with  a ^  and  ojj  the  Pedersen  and  parallel  conductivities  respectively) 
the  field  lines  may  be  regarded  as  equipotentials  and,  as  a  result,  a 
field  line  integrated  (two-dimensional)  model  is  justified  [ Volk  and 


and  Haerendel,  1971 ,  Perkins  et  al . ,  1973  ]  . 

3y  adopting  a  Cartesian  coordinate  system  (x,y,z  )  with  magnetic  field 
,  ambient  electric  field  E^y,  and  after  transforming  to  a  frame  drifting 
with  velocity  =  (cEQ/B)  x,  the  two  dimensional  model  equations  for  the 
magnetic  field  line  integrated  plasma  cloud  Pedersen  conductivity  E(x,y) 
and  the  self-consistent  plasma  cloud  electrostatic  potential  &<p(x,v) 
can  be  written 

6 1  c  "  2 

_  +  -  z  x  7<Sc?  •  VE  =  D  V  Z  (1) 

at  b  ~ 

V-  (£7oCp)  =  E  •  V£  (2) 

~o 

where  -  -  E  (x,y)  +  E  ,  E  (x,y,)  the  total  electrostatic  field,  c  is 

—  — o 

the  speed  of  light,  Vh  (d/dx,  d/dy)  and  D  is  the  cross  field  diffusion 
coefficient  [Verkins  et  al.,  1973]  given  by  2( / ft)  (ck^  T/eB)  with  T 
the  ion  and  electron  temperature  the  sum  of  electron  collision  frequen¬ 
cies  with  cloud  ions  and  ambient  neutrals,  kg  is  Boltzmann’s  constant, 
and  ft  the  electron  gyrofrequency .  For  barium  plasma  clouds,  typical 
values  of  D  lie  in  the  range  of  0.6  to  6  m^/sec  [McDonald  et  al.,  1978]  . 
All  other  symbols  retain  their  conventional  meaning.  Equation  (1)  results 
from  the  field-line  integration  (along  the  z-direction)  of  the  ion  con¬ 
tinuity  equation  while  equation  (2)  is  derived  from  current  conservation 
V- J=0. 

By  linearizing  (1)  and  (2),  i.e.,  E  =  Eq  +  <57 ,  etc.,  and  assuming  fluc¬ 
tuations  in  magnetic  field  line  integrated  Pedersen  conductivity  and  cloud 
potential  <5E,  <5  co«  exp  [i  (k^y  +  k^x)  *  k*2i  =  0,  k  L  >>1,  one 

finds  the  usual  E  x  B  growth  rate 


Y,  =  (cE  /BL)(k  /k)  -  Dk 

_k  o  y 


(3 


2  0  2  -1 

where  k  =  k  "  +  k  ,L  =  Sin  £  /dx.  For  cF  /B  =  100  m/sec,  L  =  6  km, 
x  y  o  x  o 

2 

D  =  1  m  /sec,  the  critical  wavelength  (>^=0)  “  60  m# 

Equations  (1)  and  (2)  can  be  put  in  dimensionless  form  f McDonald 
et  al .  ,  1978;  198Cl  by  normalizing  x  =  (x,y),  t,  £,  V,  <5cp  by  L0,  L0/V0, 

Z  ,  V  ,  L  E  ,  respectively,  giving 

o  o  o  o 

d_l  +  i  X?6©  -?E=R_172E 

d  t 

V  .  (E  V6  cp)  =  3E/<9y  (5; 

where  R  =  V  T,  /D.  As  a  result,  the  evolution  of  a  plasma  cloud  is  com- 
o  o 

pletely  determined  by  the  initial  cloud  configuration,  boundary  conditions 
and  the  dimensionless  number  R. 

3.  NUMERICAL  SIMULATIONS 


Equations  (1)  and  (2)  were  solved  numerically  on  a  mesh  consisting  of 
258  grid  points  in  the  x  direction  (E^  x  _B  direction)  and  102  points  in  the 
y  direction.  With  a  constant  grid  spacing  of  310  meters,  the  real  space 
dimensions  of  the  mesh  were  80  km  along  x  and  31  km  along  y.  The  magnetic 
field  line  integrated  Pedersen  conductivity  £(x,y)  in  (1)  was  advanced  in 
time  using  a  multi-dimensional  flux-corrected  variable  time  step  leapfrog- 
trapezoid  scheme  [Zalesak,  1 97 9 J  which  is  second  order  in  time  and  fourth 


order  in  space.  At  each  timestep  the  self-consistent  electrostatic 

potential  6 cp  due  to  the  ion  cloud  was  determined  from  (2)  using  a  Chebvchev 

iterative  method  [  Varga ,  1962  ;  McDonald ,  1980]  with  a  convergence  criterion 
-4 

of  10  .  Periodic  boundary  conditions  were  imposed  along  the  v  direction 

with  Neumann  conditions  (  dfd  x  =  0)  along  the  x  direction.  These  boundarv 

conditions  result  in  a  realistic  representation  of  plasma  inf low-outf low  at 

the  boundaries  in  the  E0  x  K  direction. 

The  principal  diagnostics  used  to  monitor  the  evolution  of  the  plasma 

cloud  are  the  time  history  of  the  field  line  integrated  Pedersen  conductivity 

of  the  ion  cloud  E(x,y,t),  associated  spatial  power  spectra,  and  the  plasma 

cloud  elec trostat ic  potential  6<p(x,y,t).  These  power  spectra  were  obtained 

by  first  Fourier  transforming  the  real  space  cloud  Pedersen  conductivity 

6I(x,y)  <$E(k  ,  k  )  where  6£(x,y)  =  E(x,y)-£  ,  Z  the  maximum  cloud  con- 

x  y  °2 

ductivity.  The  power  spectral  density  |6E(k  ,k)/E  |  was  then  formed  and  the 

x  '  o 

one-dimensional  power  spectra  PCk^)  and  p(k^)  were  computed  where 


P(kx)  =  /dky  |6£(kx,ky)/  T.0\ 


and 


p(ky)  =/dkx  |5E(kx,ky)  /E0| 


These  transverse  averaged  power  spectra  P(k^)  and  P(k^)  were  then  fitted 

pCeskinen  et  al . ,  1980a]  with  a  three  parameter  (spectral  strength  P  , 

spectral  index  n  ,  and  outer  scale  wavenumber  k  )  power  1  aw  of  the  form 
a  oa 


P(k  )  =  P  (1  +  k  /k  )~nc/2 
a  oa  a  oa 


where  a  =  x  or  y.  The  method  used  to  extract  the  best  fit  parameters 


5 


P  , n  ,  k  ,  is  a  nonlinear  least  squares  procedure  [Keskinen  et  al., 
oa  a  oa  L -  - 

1980a]  which  computes  P  and  n  directly  and  then  iterates  to  find  k 

J  o«  a  oa 

Initially,  the  field  line  integrated  Pedersen  conductivity  of  the 
plasma  cloud  was  taken  to  be  of  the  form 

E(x,y,t=o)&o=[M  exp(-x2/L2)  9*  0.1]  (l+e(x,v)) 

where  both  the  mean  amplitude  and  spatial  distribution  of  e(x,y)  were 
varied.  In  previous  numerical  studies  [  McDonald  et  al ■ ,  1978:  1980; 
Keskinen  et  al. ,  1980a]  e(x,y)  was  given  a  root  mean  square  value  of  3% 
and  generated  from  a  randomly  phased  Gaussian  power  spectrum.  In  this 
work  two  different  forms  for  e(x,y)  were  used.  In  the  first  case  e(x,y)  = 

A  cos  3  kpy,  a  single  monochromatic  wave  along  the  linearly  most  unstable 
y-direction  with  wavelength  2  iT/3t  =  10km  where  =  2  7r/30km  ^ . 

In  the  second  case,  r(x,y)  =  A  (l-2r(x,v))  where  r(x,y)  is  a  random  num¬ 
ber  between  0  and  1.  This  case  models  the  many  wave  initial  condition.  In 
both  cases  small  and  large  amplitude  initial  perturbations  were  modeled  by 
taking  A  equal  to  0.03  and  0.15  respectively.  Figure  1  gives  a  rough 

schematic  plot  of  the  structure  of  e(k  ,  k  )  used  in  these  simulations.  For 

x  y 

both  forms  and  amplitudes  of  e(x,y),  three  simulations  were  made  distin¬ 
guished  by  different  initial  field  line  integrated  cloud  Pedersen  conduc¬ 
tivity  scale  lengths  L  =  4,6,10  km  perpendicular  to  the  magnetic  field.  For 

all  cases  V  =  100  m/sec  and  the  cross  field  diffusion  coefficient 
o 

2 

D  =  1  m  / sec.  In  addition,  in  all  cases  M  =  1  so  that  the  maximum  field 
line  integrated  Pedersen  conductivity  of  the  cloud  to  the  background  iono¬ 
sphere  is  approximately  10.  For  the  sake  of  consistency  and  brevity,  we 
confine  ourselves  to  the  description  of  the  evolution  of  the  L  =  6  km 


barium  plasma  cloud  only. 


4.  RESULTS 


Figures  2a-2d  illustrate  the  evolution  >f  an  initially  slablike  barium 

plasma  cloud  driven  unstable  by  a  3%  monochromatic  initial  wave  by  showing 

the  real  space  isodensity  contours  of  the  field-line  integrated  Pedersen 

cloud  conductivity  in  the  plane  perpendicular  to  the  magnetic  field. 

Figure  2a  displays  the  initial  configuration.  Figure  2b  shows  the  plasma 

cloud  conductivity  at  t  =  200  sec  where  some  elongation*  steepening,  and 

jetting  to  the  frontside  has  begun.  Figure  2c  details  the  cloud  structure 

at  t  =  1000  sec  where  the  characteristic  fingers  have  formed  and  stretched 

along  the  x  _B  drift  direction.  Figure  2d  gives  the  structure  in  the 

late  time  nonlinear  regime  at  t  =  2400  sec  (y  t~  40  with  y  the  maximum 

max  max 

linear  growth  rate).  Similiar  shapes  and  morphologies  are  observed  and  for 
the  other  two  Pedersen  conductivity  gradient  scale  lengths  studied,  i.e., 

L  =  4  and  10  km,  but  on  different  time  scales. 

Figure  3  gives  representative  one-dimensional  plasma  cloud  Pedersen 
conductivity  spatial  power  spectra  computed  in  the  nonlinear  late  time 
(t  =  2400  sec)  regime  both  parallel  P(k^)  an<3  perpendicular  P(k^)  to  the 
plasma  cloud  drift  for  the  3%  monochromatic  initial  conditions  with  L  =  6  km. 


Before  significant  bifurcation  occurs,  the  perpendicular  power  spectrum 


P(ky)  is  not  described  by  a  power  law,  but  is  dominated  by  the  initially  ex¬ 
cited  mode  k  /k_  =  3  where  k_  =  (27r/30)km  *  is  the  fundamental  mode  number 
y  Fv  Fy 

in  the  y-direction.  However,  due  to  the  steep  edges  of  the  plasma  cloud 


striations  in  the  parallel  x-direction  the  parallel  power  spectra  P(k  ) 
do  conform  to  a  power  law  with  spectral  index  of  approximately  2.  At 


late  times  in  the  nonlinear  regime,  however,  the  perpendicular  P(k^)  power 

spectra  also  tend  toward  a  power  law. 

Figure  4  shows  the  time  histories  of  n  and  n  for  the  small  amplitude 

xv 

monochromatic  initial  conditions  with  L  =  6  km.  Similar  results  were  found 
for  L  =  4,  10  km. 


Figure  5  gives  a  sample  contour  plot  of  the  unaveraged  two-dimensional 
spatial  power  spectra  of  the  plasma  cloud  Pedersen  conductivity  at  t  =  2400 
sec  using  3%  monochromatic  initial  conditions  with  L  =  6  km.  As  can  be 
seen,  the  two-dimensional  spatial  power  spectra  is  anisotropic  with  most  of 
the  spectral  power  concentrated  along  and  near  the  linearly  most  unstable 
v-direction.  This  anisotropy  in  (k^k^)  space  is  consistent  with  the  late 
time  distribution  in  (x,y)  space  of  the  field  line  integrated  cloud  Pedersen 
conductivity  which  show  striations  elongated  muGh  more  in  the  x-direction 
OEq  x  .B  drift  direction).  Similar  anisotropic  spatial  power  spectra  in  the 
late  time  nonlinear  regime  are  also  found  using  the  large  amplitude  (15%) 
monochromatic  initial  conditions. 

Figures  6a-6d  give  the  real  space  distribution  of  field  line  integrated 
cloud  Pedersen  conductivity  for  the  large  amplitude  15%  monochromatic  initial 
conditions  for  L  =  6  km.  Figure  6b  describes  the  conductivity  profile  also 
at  t  =  200  sec  and  shows  a  similar  configuration  to  the  3%  case  in  Figure  2b 
with  the  exception  that  more  elongation  and  jetting  to  the  frontside  has 
taken  place  in  the  larger  amplitude  case.  For  the  monochromatic  initial  con¬ 
ditions  used  the  evolution  of  the  plasma  cloud  with  the  larger  amplitude  per¬ 
turbation  was  found  to  proceed  on  a  faster  time  scale  than  the  plasma  cloud 
seeded  with  a  smaller  amplitude  perturbation.  Figure  6c  displays  the  cloud 
contours  at  t  -  1000  sec  while  Figure  6d  details  the  plasma  cloud  field  line 
integrated  Pedersen  conductivity  at  t  =  2400  sec  and  is  not  dissimilar  from 
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the  plasma  cloud  configuration  using  the  3%  monochromatic  initial  conditions 
(cf.  Fig.  2d).  Similar  morphologies  are  also  observed  for  the  olasma  clouds 
with  conductivity  gradient  scale  lengths  L=4,6,10  km,  but  on  different  time 
scales  with  the  L=4  km  case  proceeding  faster  and  the  L=10  km  case  slower. 

Figure  7  gives  sample  one-dimensional  parallel  P(k^)  and  perpendicular 
P(k^)  spatial  power  spectra  of  plasma  cloud  conductivity  computed  in  the 
nonlinear  late  time  (t  =  2400  sec)  regime  for  15%  monochroma t ic  initial  con¬ 
ditions  with  L  =  6  km.  Note  the  similarity  with  Figure  3. 

Figure  8  displays  the  time  dependence  of  the  best  fit  spectral  indices 

n  and  n  for  the  15%  monochromatic  initial  perturbation  with  L  =  6  km. 
x  y 

As  can  be  noted  the  spectral  indices  n^  and  achieve  their  quasi-steady 
state  values  on  a  slightly  faster  time  scale  (on  the  order  of  100-200  sec) 
for  the  larger  amplitude  monochromatic  perturbation  a?  opposed  to  the  smaller 
amplitude  monochromatic  initial  conditions  (see  Figure  4).  The  value  for  the 
spectral  index  in  the  ExB  direction  approaches  n^  -  2  while  n^  -  2-3. 

Similar  spectral  indices  are  also  found  for  the  other  two  cases  studied,  i.e., 
L  =  4,  10  km,  but  again  on  different  time  scales. 

Figures  9a-9d  portray  the  evolution  of  the  plasma  cloud  using  purely 
random  initial  conditions  with  maximum  amplitude  of  3%  for  L  =  6  km.  In 
these  simulations  the  initial  evolution  of  a  test  wave  is  influenced  not  onlv 
by  the  ambient  plasma  cloud,  but  also  by  a  many  wave  background.  Figure  9b 
displays  the  cloud  at  t  =  200  sec  and  shows  the  initial  random  perturbations 
developing  on  the  backside.  Figure  9c  illustrates  the  field  line  integrated 
Pedersen  conductivity  contours  at  t  =  1000  sec  with  striation  and  elongation 
evident.  Figure  9d  shows  the  cloud  configuration  at  t  =  2400  sec  and  is 
similar  to  the  late  time  configurations  using  the  monochromatic  initial  con¬ 
ditions  as  shown  in  Figure  2d  and  6d.  The  evolution  of  the  clouds  with 
gradient  scale  lengths  L  =  4,10  km  under  3%  random  initial  conditions  is 


similar  to  the  L  =  6  km  case,  but  on  different  time  scales,  i.e.,  the  L  =  4 
km  cloud  structuring  faster  in  absolute  time  and  the  L  =  10  km  cloud 
developing  slower. 

Figure  10  shows  sample  one-dimensional  power  spectra  of  plasma  cloud 

conductivity  with  L  =  6  km  both  parallel  P(k  )  and  perpendicular  P(k  )  to 

x  y 

the  ExB  drift  (x-direction)  at  t  =  2400  sec  for  the  purely  random  3%  initial 

conditions.  The  time  histories  of  the  best  fit  spectral  indices  n  and  n 

xv 

both  in  the  parallel  P(k  )  and  perpendicular  P(k  )  directions  with  L  =  6 

x  y 

km  for  the  small  amplitude  (3%)  random  initial  conditions  are  displayed  in 

Figure  11.  After  initial  transients,  the  spectral  index  in  the  x  jl 

direction  becomes  n^  -  2  with  ~  2-3.  In  comparing  the  time  evolution  of 

n  and  n  for  the  3%  monochromatic  and  random  initial  conditions  (see  Figures 
x  y 

4  and  11)  one  notes  that  the  spectral  indices  reach  their  quasi-steady  state 
values  faster  for  the  many  wave  random  initial  conditions.  This  is  partic¬ 
ularly  true  for  the  spectral  index  n^.  The  same  values  for  the  spectral 
indices  and  n^  are  also  observed  using  the  gradient  scale  lengths 
L  =  4,10  km  starting  from  the  3%  random  initial  conditions.  These  spectral 
indices  are  in  agreement  with  both  experimental  values  [Baker  and  Ulwick, 

1978;  Kelley  et  al . ,  1979]  and  previous  one  level  [McDonald  et  al . , 

1980;  Keskinen  et  al . ,  1980]  and  two-level  [Scannapieco  et  al . ,  1976] 
numerical  simulations  of  ionospheric  barium  clouds  using  different  initial 
conditions . 

The  turnover  or  outer  scale  size  in  the  one-dimensional  perpendicular 
power  spectra  P(k^)  of  the  plasma  cloud  conductivity  with  L  =  6  km  using  the 
small  amplitude  (3%)  random  initial  perturbations  is  plotted  in  Figure  12  as 
a  function  of  time.  In  the  early  nonlinear  regime,  the  perpendicular  outer 
scale  size  2tt / k^^  approximates  the  initial  parallel  Pedersen  conductivity 

/ 
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gradient  scale  length  L  of  the  cloud  in  agreement  with  previous  numerical 
simulations  [Keskinen  et  al . ,  1980a]  using  Gaussian  initial  conditions 

with  mean  amplitude  of  3%. 

Figure  13  shows  a  representative  sample  of  the  unaveraged  two-dimensional 
spatial  power  spectra  of  the  cloud  Pedersen  conductivitv  of  t  =*  2400  sec 
using  the  3%  random  initial  conditions  with  L  =  6  km.  The  power  spectra 
is  anisotropic  with  most  of  the  power  located  along  and  near  the  linearly 
most  unstable  y-direction.  Note  the  similarities  with  Figure  5  which 
evolved  from  monochromatic  initial  conditions.  Similar  anisotropy  is  found 
for  L  =  4,10  km. 

Figures  14a-14d  display  the  evolution  of  the  barium  plasma  cloud  using 
the  15%  purely  random  initial  conditions  for  the  L  =  6  km  scale  length.  Com¬ 
paring  the  plasma  cloud  structure  in  Figure  14b  at  t  =  200  sec  with  Figure 
9b  (small  amplitude  random  initial  perturbation  at  t  =  200  sec)  one  again 
notes  the  increased  elongation  and  penetration  to  the  frontside  in  the  larger 
initial  amplitude  case.  Figure  14c  shows  the  striations  at  t  =  1000  sec 
while  Figure  14d  illustrates  further  bifurcation  at  t  =  2400  sec.  Many 
similarities  can  be  seen  in  comparing  Figure  14d  and  Figure  9d. 

Figure  15  shows  sample  one-dimensional  power  spectra  of  the  L  =  6  km 
plasma  cloud  Pedersen  conductivity  both  parallel  P(k  )  and  perpendicular 
P(k  )  to  the  ExB  drift  (x-direction)  at  t  =  2400  sec  for  the  purely  random 

y  -  ’ 

15%  initial  conditions.  The  time  histories  of  the  best  fit  spectral  indices 

n  and  n  both  in  the  parallel  P(k  )  and  perpendicular  P(k  )  directions  with 
x  y  x  v 

L  =  6  km  for  the  large  amplitude  (15%)  random  initial  conditions  are  display¬ 
ed  in  Figure  16.  The  spectral  index  in  the  E  x  B  direction  (x-direction) 

— o  — 

approaches  -  2  with  the  spectral  index  perpendicular  to  the  E^  x  B 
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drift  n ^  -  2-3.  By  comparing  Figure  16  (15%  random  initial  conditions)  with 

Figure  11  (3%  random  initial  conditions)  one  again  notes  that  the  spectral 

indices  n  and  n  approach  their  quasi-steadv  state  values  on  a  sllghtlv 
x  y 

faster  time  scale  for  the  larger  amplitude  initial  perturbations.  In  ad¬ 
dition,  in  comparing  Figure  16  (15%  random  initial  conditions)  with  Figure  8 
(15%  monochromat ic  initial  conditions)  it  can  be  said  that  the  spectral  in¬ 
dices  n  and  n  develop  on  a  faster  time  scale  for  the  mnnv  wave  initial  per- 
xv 

turbations  as  opposed  to  the  single  wave  monochroma t i c  initial  conditions. 

The  turnover  or  outer  scale  size  In  the  one  dimensional  perpendicular 

power  spectra  P(k^)  for  the  L  =  6  km  plasma  cloud  seeded  with  the  large 

amplitude  (15%)  random  initial  perturbations  corresponding  to  Figure  14  is 

plotted  in  Figure  17  as  a  function  of  time.  In  the  nonlinear  regime  (for 

the  times  the  simulations  were  run)  the  perpendicular  outer  scale  size 

2n/k  approximates  the  initial  parallel  (to  E  x  B  motion)  Pedersen 
oy  — o 

conductivity  gradient  scale  length  L  of  the  cloud  in  agreement  with  the  3% 
random  initial  conditions  (see  Figure  12).  Again,  this  scaling  is  achieved 
on  a  slightly  faster  time  scale  than  with  the  smaller  amplitude  3%  random 
initial  perturbations. 

5.  SUMMARY  OF  RESULTS 

We  have  studied  the  effects  of  different  initial  conditions  on  the 
evolution  of  large  F  region  ionospheric  plasma  clouds  driven  unstable  bv  the 
ExB  gradient  drift  instability .  This  has  been  accomplished  bv  the  numerical 
solution  of  the  fundamental  one  level  (one  for  the  plasma  cloud  and  neglect 
of  cloud -background  ionosphere  interaction)  two-dimensional  magnetic  field 
line  integrated  plasma  cloud  fluid  equations  using  several  different  sets  of 
initial  conditions.  In  these  numerical  simulations  both  large  and  small 


purely  random  and  monochromatic  initial  conditions  have  been  used  ?<>  m 
initially  slab  like  plasma  cloud  models  with  different  magnetic  tit  Id  !i 
tegrated  Pedersen  conductivity  gradient  scale  lengths  I  =  4,^,1 o  km  *»  r 
dicular  to  the  magnetic  Held.  Independent  of  the  initial  condition 
ied,  all  models,  in  the  nonlinear  late  time  regime  become  similar  hot 
real  (x,v)  space  and  in  Fourier  (k  ,  k^  1  space  where  x  and  \  det  itn  tne 
perpendicular  to  the  direction  of  the  magnetic  field.  In  the  non  I  i  :n  if 
gime,  an  outer  scale  turmwer  is  observed  in  the  spatial  power  spectr.i 
puted  in  a  direction  (v)  perpendicular  to  the  l\xK  drift  ( v )  of'  the  :*1  i 
cloud.  Under  the  different  initial  conditions  studied,  tne  tvo-  i 

spatial  power  spectra  in  the  nonlinear  reel  me-  in  rh«*  plane  riu-n.li.  i! 
the  magnetic  field  is  found  to  be  anisotropic-  with  most  ,*!  the  sr«-»  t  r  ■ 
concentrated  in  the  linear  lv  most  unstable  v-dirrt't  jim.  Furthermore, 

independent  of  the  initial  conditions  used  in  this  report  the  one-dimer 

-nx 

x  power  spectra  is  «  k  with  n  2  for  J  "  between  1  and  Hu  kr, 

XX  X 

the  v  power  spectra  «  k  v  with  nv  J- 1  for  J  -  ;k  between  1  and  10k«' 

v  v 

the  random  initial  perturbation  case  achieves  its  quasi-final  power  law 
spectrum  description  on  a  slightly  faster  time  scale-  than  the  sine  It  mo 
chromatic  wave  initial  perturbation  case.  Within  the  context  of  initi  \ 
perturbation  amplitude,  for  a  given  case,  the  larger  amplitude  nertur ■  i 
attains  the  quasi-final  state  more  rapidlv.  These  numerical  result^  « 
consistent  with  recent  experimental  Qftaker  and  Ulwjck,  1978;  Kell*  if 
1976;  Chat urved  1  and  Ossakow ,  1979;  Kesk_inen  et  al  .  ,  1980^  studies  • 
Ionospheric  plasma  clouds. 

Future  analytical  and  numerical  studies  of  the  evolution  of  ionoev 
plasma  clouds  are  planned  which  include  the  addition  of  inertial  effect 
the  coupling  to  other  ionospheric  levels. 
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y),  Curve  A  refers  to  previous  simulations  [McDonald  et  al.,  1978; 
has  a  Gaussian  dependence  on  k  =  (kx2  +  ky2)1/2.  Curve  B  repre¬ 
curve  C  is  approximately  flat  in  k. 


A 


1000  SEC 


(c)  (d) 

Fig.  2  —  Real  space  isodensity  contour  plots  of  E(x,y)/E0  for  L  =  6  km 
using  3%  monochromatic  initial  conditions  at  (a)  t  =  0  sec,  (b)  t  =  200  sec, 
(c)  t  =  1000  sec,  and  (d)  t  =  2400  sec.  Eight  contours  are  plotted  in  equal 
increments  from  0.1  to  1  with  every  other  contour  represented  by  a 
dashed  line.  The  x  axis  (y  axis)  denotes  the  E0  X  B(E0)  direction. 
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Fig.  5  —  Contour  plot  of  log10  l6S(kx,ky  )/£0  I2  at  t  =  2400  sec  for  L  =  6  km  using  the 
3%  monochromatic  initial  conditions.  Six  contours  are  plotted  which  correspond  to 
decades  of  decreasing  spectral  power.  The  values  of  kpx  =  27T/80  km'1  with  kFv  = 
27T/31  km"1. 


300  600  900  1200  1500  1800  2100  2400  t(sec) 

Fig.  8  —  Time  history  of  best  fit  spectral  indices  nx  and  ny  for  L  =  6  km  using  the 
15%  monochromatic  initial  conditions.  The  spectral  index  ny  is  not  plotted  before 
approximately  1200  sec  since  P(ky )  is  not  a  power  law. 
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Fig.  10  —One  dimensional  (a)  x  power  spectra  P(kx)  and  (b)  y  power  spectra  P(ky)  at  t  =  2400  sec 
for  L  *  6  km  using  the  3%  random  initial  conditions  with  n_  =  2.0,  2x/k„_  =  32  km,  and  n„  =  2.4, 


koyL/47r  for  4,  6,  10  km  using  the  3%  ra 
ns.  7(L)  is  maximum  linear  growth  rate. 


(c)  (d) 

Fig.  14  —  Real  space  isodensity  contour  plots  of  l(x,y  )/I€)  for  L  =  6  km  using  the 
15%  random  initial  conditions  at  (a)  t  *  0  sec,  (b>  t  *  200  sec,  (c)  t  *  1000  sec,  and 
(d)  *  2400  sec 
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Fig.  15  —  One  dimensional  (a)  x  power  spectra  P(kx)  and  (b)  y  power  P(ky)  at  t  *  2400  sec 
for  L  *  6  km  using  the  15%  random  initial  conditions  with  nx  *  2.0,  2ir/kox  *  30  km  and 
n„  -  2.3,  2jr/k_v  ■  3.9  km 
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Fig.  17  —  Time  history  for  koy  L/4n  for  L  =  4,  6, 10  km  using  the  15%  random 
initial  conditions.  >(L)  is  maximum  linear  growth  rate. 
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OICY  attn  DlNER  8AR?_£T^  *?1 

OICY  ATT\  R5 


DEFENSE  COMMUNICATION  ENGINEER  CENTER 
1860  wIEHLE  AVENUE 
RESTON,  VA.  22090 

OICY  attn  CODE  fR 820 

OICY  ATTN  COOE  R410  JAMES  W.  MCLEAN 

OICY  attn  COOE  R720  J.  WORThI NGrON 

DIRECTOR 

DEFENSE  COMMJNI CAT IONS  AGENCY 
WASHINGTON,  D.C.  20305 

(ADR  CNWDI:  attn  COOE  240  FOR) 
OICY  attn  COOE  1016 


COMMANDANT 

NATO  SCHOOL  (SHAPE. 

APO  NEW  YORK  091/2 

31CY  attn  j.5.  DOCUMENTS  CF-1CER 

UNDER  SECY  OF  DEF  FOR  RSCH  i  EnGRG 
DEPARTMENT  OF  DEFENSE 
WASHINGTON,  D.C.  20301 

OICY  attn  STRATEGIC  i  SPACE  SYSTE 

WWMCCS  SYSTEM  ENGINEERING  CRG 
WASHINGTON,  D.C.  20305 

OICY  attn  R.  CRAWFORD 


DEFENSE  TECHNICAL  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA.  22314 

(12  COPIES  IF  OPEN  PUBLICATION,  OTHERWISE  2 
12CY  ATTN  TC 


COMMAM)ER/DIRECtOR 
ATMOSPHERIC  SCIENCES  MORATORY 
CQP{£S)  J.S.  ARMY  ELECTRONICS  COMMAND 

white  SANCS  MISSILE  RANGE,  nm  3800 2 
OICY  attn  DE^AS-cO  F.  NI.E5 


DIRECTOR 

DEFENSE  INTELLIGENCE  AGENCY 
WASHINGTON,  O.C.  20301 


OICY 

ATTN 

OT-18 

OICY 

attn 

D8-4C 

E. 

0,cARRELL 

OICY 

ATTN 

DIAAP 

A. 

WISE 

OICY 

attn 

DIA5T- 

5 

OICY 

ATTN 

0T-19Z 

R, 

MORTQN 

OICY 

ATTN  HQ-TR 

J. 

STEwART 

OICY 

ATTN 

w.  wITtig 

DC-7D 

DIRECTOR 

DEFENSE  NUCLEAR  AGENCY 
WASHINGTON,  D.C.  20305 
OICY  ATTN  STVL 
04CT  ATTN  TITL 
OICY  Attn  DOST 
05CY  ATTN  RAAE 

COMMANOER 
field  COMMAND 
DEFENSE  NUCLEAR  AGENCY 
KIRTlAND  AFB,  NM  87115 
OICY  ATTN  FCPR 

director 

INTERSERVICE  NUCLEAR  WEAPONS  SCHOOL 
KIRTLAND  AF0,  ^  87115 

OICY  ATTn  DOCUI-eNT  CONTROL 


DIRECTOR 

3MD  ADVANCED  TECH  CtR 
HUNTSVILj-E  OFFICE 
P.  O.  BOX  1500 
HUNTSVILLE,  Al  35807 

OICY  attn  ATC-T  MELVIN  CAPPS 
OICY  attn  ATC-G  W.  DAVIES 
OICY  attn  ATC-R  DON  RUSS 

program  manager 

3MD  PROGRAM  OFFICE 
5001  E ISENhOWER  AVENUE 
ALEXANDRIA,  VA  22553 

OICY  ATTN  DACS-3M"  v-E-* 

CH/Er  C-E  SERVICES  DIVISION 
U.S.  ARMY  COMMUNICATIONS  CMC 
PENTAGON  RM  13269 
WASHINGTON,  D.C.  <0510 

olcy  attn  :-e-ser\ i :es  dw;5;:n 

commander 

FRAOCOM  Technical  Support  aCtn;'v 
DEPARTMENT  OF  t-,£  army 
FORT  mONMOuTh,  n.j.  07703 

OICY  attn  DRSEl-nl-RD  BE'VwE" 
OICY  attn  DRSEL-PL-ENV  h,  9oMC£ 
OiCY  AT tn  j .  F  .  QulGuEy 


COMMANDER 

HARRf  DIAMOND  .A0ORA  TORIES 
Department  jf  Th£  army 
28QQ  ^OwOER  MIll  ROAD 
AC£lPh{,  md  20735 

(CNwDl-INI^R  ENVELOPE :  ATTn: 
0 ICY  Atjn  ^hO-TI  M.  wEINER 
0  ICY  ATTn  DELmD-RB  R.  w  III l  AMS 
0 ICY  ATTN  DELHD-NP  r.  WlMENlTZ 
0  ICY  ATTN  DElmD-NP  C.  MOAZED 


DELMD-RBH) 


commander 

J.S.  ARMY  CQMM-ELEC  ENGRG  INStAl  AGY 
ct.  -uACMuCA,  A Z  35613 


0  ICY  ArTN  CCC-cMEO  GEORGE  LANE 


QIRECTOR 

NAVAL  RESEARCH  ^A0ORAtO«Y 
WASHINGTON,  D.C ,  20375 


OICY 

ArTN  COOE 

4700 

T.  P.  COFfEr  (25  CYS  IF  JN 

OICY 

ATTN  code 

4701 

JACK  D.  BROWN 

OICY 

attn  COOE 

4780 

3 RANCH  HEAD  (150  CYS  IF  UN, 

OICY 

attn  COOE 

7500 

HQ  COMM  DIR  BRUCE  wAlD 

oicr 

ArrN  COOE 

7550 

J.  DAVIS 

OICY 

ArTN  CODE 

7580 

OICY 

ArTN  COOE 

7551 

oicy 

A^n  COOE 

7555 

OICY 

ATTN  CODE 

4730 

E.  MCLEAN 

COMMANDER 

'NAVAL  SEA  5YST£.MS  COMMAND 
WASHINGTON,  D.C.  20562 

01CY  ATTN  CApT  R.  PI^JN 


COMMANDER 

-.3.  ARMY  -OREIGN  SCIENCE  &  TECH  CTR 
220  ?Tr.  5tREEt,  NE 
^ARLOrT£Sv!-^E,  VA  22901 
0 ICY  attn  ORXST-50 
01CY  ATT^  R.  ^ONES 


COMMANDER 

-j.3.  ARMY  MATER 1EL  DEV  S  READINESS  CMD 
500 1  EISENhOER  AVENUE 
ALEXANDRIA,  VA  22335 

3  ICY  ATtn  ORCwX  J.  A.  9ENOER 


commander 

. . c .  army  nuclear  ano  chemical  agency 
7500  3AOQUIOC  ROAD 
sl'jg  :o^ 

SPRINGFIELD.  vA  22150 
01CY  ArTN  LIBRARY 

DIRECTOR 

j.S.  ARMY  3AlLIStIC  RESEARCH  uA8S 
ABERDEEN  MOVING  GROUND,  MD  21005 

01CY  attn  tech  lIB  EDWARD  BA ICY 

COMMANDER 

,  .  5  .  ARMY  SATO  CM  AGENCY 
=  ~.  HONMOuTH,  nj  07705 

5icr  attn  document  control 

COMMANDER 

j.S.  army  MISSILE  INTELLIGENCE  AGENCY 
REDSTONE  arsenal,  AL  35809 
0 ICY  ATTN  JIM  GAMBLE 

director 

..5.  ARMY  rRADOC  SYST£MS  ANALYSIS  ACTIVlTY 
WHITE  SANDS  "ISSILE  RANGE,  nm  38002 
OICy  ATTN  ATAA-SA 
J 1 C Y  ATTN  rCC/P.  PAYAN  uR. 

OiCY  AT^N  ATAA-TAC  LTC  d.  HESSE 


COMMANDER 

naval  space  surveillance  system 

DAHLGREN,  VA  22448 

01CY  ArTN  CAPT  j.  H.  BURTON 

officer- in-charge 
naval  surface  weapons  center 
white  oax,  silver  spring,  md  20910 
0 ICY  AttN  CODE  F3I 

OIRECTOR 

STRATEGtC  SYSTEMS  ^ROdEC  OFFICE 
DEPARTMENT  QF  tH£  navy 
wASHtNGTON,  O.C.  20376 
01CY  ATtn  NSP-2141 
01CY  Attnj  nSSP-2722  FRED  WIMBERLY 

NAVAL  SPACE  SY$TEM  AdV^Y 
P.  O.  BOX  92960 
WORLDWAY  PQStAl  CENTER 
lOS  ANGELES,  CALIF.  90009 

01CY  ATTN  A.  B.  HAZZARD 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CENtFR 
DAMLGREN  LABORATORY 
DAhLGREN,  VA  22448 

01CY  attn  CODE  OF -14  R.  9UTLER 

COMMAND l NG  OFFICER 

NAVY  SPACE  SYSTEMS  ACTtVlTY 

p.O.  SOX  92960 

WORLDWAY  postal  CENTER 

wCS  ANGELES,  CA.  900Q9 
01CY  ATtn  CODE  52 

OFFICE  OF  NAVAl  RESEARCH 
ARLINGTON,  VA  22217 

01CY  AttN  CODE  465 
01CY  Attn  COOE  461 
01CY  Attn  CODE  402 
01CY  A"n  CODE  420 
01CY  ATtn  CODE  421 


;OM MANOER 

naval  electronic  sy$tems  ccmmand 
WASHINGTON,  D.C,  20360 

0  ICY  ATTN  nAVAlEX  034  T.  HUGHES 
OICY  ATTN  ^ME  117 
01CY  attn  PME  1X7  — t 
OICY  ArTfv|  COOE  5011 

CCmANOlNG  OFFICER 
naval  intelligence  support  crR 
4301  SUrLAND  ROAD.  3LQG.  5 
wAShIn&tON,  D.C.  20590 

OiCY  attn  dubbin  stic  12 

oicy  at^n  n isc -50 

OICY  ArTN  CODE  5404  J.  GAU£T 


_OMMAnOER 

naval  OCEAN  SYSTEMS  CENrER 
SAN  OlEGO,  CA  92L52 

0JCY  A^N  CODE  532  w.  MOlER 
OICY  attn  COOE  3250  C.  8AGGETT 
jICY  Attn  CODE  81  R .  EASTMAN 


COLANDER 

AEROSPACE  DEFENSE  COMMAND /DC 

department  OF  thE  AIR  corce 

ENT  AF3f  CQ  g0gi2 

OICY  ATTN  DC  MR.  long 

commander 

AEROSPACE  DEFENSE  COMMAND/ XPD 
department  OF  tH£  A(R  FORCE 
ENT  AF8 ,  CO  80912 

OICY  ATtN  xPDQQ 

02CY  attn  xp 

AIR  FORCE  GEOPHYSICS  LABORAtO«Y 
HANSCOM  AF9,  MA  01731 

oicy  attn  opr  harold  Gardner 

OICY  A*tn  OPR-l  JAMES  C.  ULWlCK 

OICY  Attn  l<b  XENNEth  S.  w.  CHAMPION 

OICY  ATtn  opr  alva  T.  StAU 

OICY  A^Tn  PHP  JULES  AARONS 

OICY  attn  PHD  JURGEN  BuChAu 

OICY  ATTN  PHD  JOHN  P.  HJLLEN 


F  CLASS) 
F  CLASS) 
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A F  WEAPONS  LAB ORA r OR Y 
KIRTLAND  AFB,  NM  87117 
0  ICY  ATTN  SUL 
Q1CY  ATTN  CA 
0  ICY  ATTN  DYC 
OICY  ATTN  DYC 
01CY  ATtn  DYT 
0  ICY  ATTN  DES 
Q1CY  ATTN  OYC 


arthur  m.  guenther 

CAPT  J.  BARRY 
JOHN  M.  <AMM 
CAPT  MARK.  A.  PRY 
MAj  GARY  G AMONG 
J.  JAW l 


AFTAC 

PATRI<X  APB,  PL  32925 

01CY  ATTN  TF/MAJ  WILEY 
01CY  ATTN  TN 

air  force  Wright  Aeronautical  Laboratories 

WRIGHT -PATTERSON  AFB,  OH  45433 
OICY  ATTN  AAD  WADE  HUNT 
01CY  ATTN  AAD  ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 

RESEARCH,  DEVELOPMENT,  J  ACQ 

DEPARTMENT  Of  THE  AIR  FORCE 

WASHINGTON,  D.C.  2033Q 
OICY  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  OIVISION/XR 

DEPARTMENT  OF  THE  AIR  FORCE 

HAMSCOM  AFB,  MA  01731 

01CY  ATTN  XR  J.  OEAS 


COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

OICY  at-^n  EEP  A.  lorentzen 

DEPARTMENT  OF  ENERGY 

ALBUQUERQUE  QPERAtIOnS  OFFICE 
P.  0.  80X  5400 
ALBUQUERQUE,  NM  87115 

OICY  Attn  DOC  CON  FOR  D.  ShERwOOO 

DEPARTl^ST  OF  ENERGY 
L18RARY  ROOM  G-042 
WASHlNGTON,  D.C.  20545 

OICY  AT^n  DOC  CON  FOR  A.  LABOWIT^ 

EG6G,  INC. 

LOS  ALAMOS  DIVISION 

P.  0.  BOX  809 

LOS  ALAMOS,  NM  85544 

OICY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.  0.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 

OICY  ATTn  DOC  CON  FOR  L-389  R.  OTT 

OICY  ATTN  DOC  CON  FOR  L-3 1  R.  HAGER 

OICY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/YSEA 

DEPARTMENT  of  the  AIR  FORCE 

HANSCOM  AFB,  MA  01731 
OICY  ATTN  YS£A 

HEADQUARTERS 

electronic  systems  division/dc 
department  of  the  air  force 

HANSCOM  AFB,  MA  01731 

OICY  ATTN  DCKC  MAj  J.C.  CLARK 

~qmfwnder 

foreign  technology  division,  AFSC 

WRIGHT -PATTERSON  AFB,  OH  45433 
OICY  ATTN  NJCD  LIBRARY 
OICY  ATTN  STOP  B.  BALLARD 

commander 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 

GRIFF ISS  AFB,  NY  13441 

OICY  ATTN  DOC  L IBRARY/TSlO 
OICY  ATTN  OCSE  V.  COYNE 

SAMSO/SZ 

POST  OFFICE  BOX  92960 

WORLDWAY  POSTAL  CENTER 

los  Angeles,  ca  90009 

(SPACE  DEFENSE  SYSTEMS) 


OICY 

ATTN 

SZj 

STRATEGIC 

AIR  COMMAND / XPF S 

OFFUTT  AFB 

,  NB 

68113 

OICY 

Attn 

XPFS  MAJ  9.  StEPmAN 

OICY 

ATtn 

ADWATE  MAJ  BRUCE  BAUER 

OICY 

attn 

NRT 

OICY 

AtTn 

DOK  CHIEF  SCIENTIST 

LOS  ALAMOS  SClENT IF IC  LA80RAT0RY 

p.  0.  BOX  1663 

LOS  ALAMOS,  NM  87545 

OICY  AT^N  DOC  CON  FOR  J.  WOLCOTT 

OICY  AT1rN  DOC  CON  FOR  R.  F.  TASCHEK 

OICY  AT^n  DOC  CCJN  FOR  £.  JONES 

OICY  A-H-n  DOC  CON  FOR  J.  MAlIK 

OICY  ATTN  DOC  CON  FOR  R.  JEFFRIES 

OICY  AT^n  DOC  CON  FOR  J.  ZlNN 

OICY  AT^n  DOC  CON  FOR  P.  KEATON 

OICY  A'rTN  DOC  CON  FOR  D.  wESTERVELT 

SAND  I A  LABORATORIES 
P.  0.  BOX  5800 
ALBUQUERQUE,  NM  87115 

OICY  A^N  DOC  CON  cOR  J.  MARGIN 

OICY  ATTN  DOC  CON  FOR  W.  BROWN 

OICY  ATTN  DOC  CON  FOR  A.  THOR N8 ROUGH 

OICY  AT^n  DOC  CON  FOR  T.  WRIGHT 

OICY  ATTN  DOC  CON  FOR  D.  DAHLGREN 

OICY  AT^N  DOC  CON  FOR  5141 

OICY  AT^n  DOC  CON  FOR  SPACE  PROUECT  01 V 

SAND I A  LABORATORIES 
LIVERMORE  <_ABORATORY 
P.  0.  BOX  969 
LIVERMORE,  CA  94550 

OICY  Attn  DOC  CON  FOR  B.  NURPMEY 

OICY  An'N  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MI L l TARY  APPLICATION 
department  of  energy 
WASHINGTON,  D.C.  20545 

OICY  ATtn  DOC  CON  FOR  D.  GALE 


other  government 


SAMSO/SK 
P.  0.  BOX  92960 
WORLOWAY  PQStAl  CENTER 
LOS  ANGELES,  CA  90009 

OICY  ATtn  SKA  (SPACE  COMM  SYSTEMS)  M. 


central  intelligence  agency 

attn  RD/S l,  RM  5G48,  HQ  BLDG 
WASHlNGTON,  D.C.  20505 

CLAVlN  OICY  AT^n  OSI/PSID  RM  5F  19 


SAMSO/W 

NORTQN  AFB,  CA  92409 

(MINUTEMAN) 

OICY  ATTn  LTC  KENNEDT 


department  of  cov-erce 
national  bureau  of  stam)ard$ 

wAShIn&tON,  D.C.  20234 

(ALL  CORRES:  ATtn  SEC  OFFICER  FOR) 
OICY  ATtn  ft.  MOORE 
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lNSrruTE  FOR  T£l£COM  SCIENCES 
NAtIONAl  TELECQMMUN 1 CAr I QN$  &  INFO  ADMIN 
BOULDER,  CO  30503 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 
OICY  ATTN  W.  JTLAlT 
01CY  ATTN  0.  CROMBtE 
0 ICY  AttN  L.  berry 

NATIONAL  OCEANIC  Sr  ATMOSPHg.RIC  ADM  IN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
department  j f  COMMERCE 
BOULDER,  CO  30302 

a  icy  a^n  r.  grubs 

31 CY  ATT,N  AERONCMY  - AB  G.  REID 


department  of  defense  contractors 


AEROSPACE  CORPORATION 
0,  BOX  92957 
-OS  ANGELES,  CA  90009 

01CY  ATtn  I.  GARFUNKEL 
01CY  ATTN  T.  SAlMI 
01CY  ATTN  V.  JOSEPHSON 

aicv  at^n  s.  sower 

01CY  ATTN  N.  STQCXWELL 
01CY  AT^N  D.  OLSEN 


01CY  AY^n  SMFA  for  ^ww 

analytical  SYSTEMS  ENGINEERING  corp 
5  OLD  CONCORD  ROAD 
BURLINGTON,  HA  01803 

0 ICY  ATTN  RADIO  SCIENCES 


BERKELEY  RESEARCH  ASSOClATES,  INC. 

P.  0.  BOX  983 
3ERKELEY,  CA  94701 

QLCY  ATTN  uCRKMAN 

BOEING  COK’ANY,  THE 
P.  0.  BOX  3707 
SEATTLE,  wA  98124 

01CY  At^n  G.  KEISTER 
QLCY  ATTN  D.  MURRAY 
01CY  AT’N  G.  MALL 
OICY  AT^N  J.  <E!*4EY 

CALIFORNIA  Ar  SAN  DIEGO,  UN IV  OF 

P.O.Box  6049 

San  Diego,  CA  92106 

BROWN  ENGINEERING  COMPANY,  INC. 
CUMMINGS  RESEARCH  3 ARK 
HUNTSVILLE,  Al  35807 

OiCY  A^n  ROMEO  A.  DEL t SERI  S 

CHARLES  StARK  DRAPER  LABORATORY,  INC, 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 

OICY  Attn  D.  3.  COX 
OICY  AT^N  J.  p.  GILMORE 


COMSAT  LABORATORIES 
Lt nth (CUM  ROAO 
CLARKSBURG,  MO  20734 

3 ICY  Ar-N  j.  ~ro£ 

CORNELL  UNIVERSITY 

department  .jf  electrical  engineering 
ITHACA,  NY  14850 

OiCY  A^n  0.  T.  FARLEY  JR 


ELECTROSPACE  SYSTEMS,  INC. 

SOX  1359 

RICHARDSON,  T K  75080 

OICY  At^n  h.  uOGSrON 

OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESL  INC. 

495  JAVA  DRIVE 

Sunnyvale,  ca  94086 

OICY  Attn  j,  ROBERTS 
OICY  Attn  JAMES  MARSHALL. 

OICY  AT  tv  c.  w.  PR£TTIE 

FORD  AEROSPACE  S  COmjN !  CAT I  ONS  CORP 
3939  FABIAN  *AY 
pALO  AlyO,  CA  94305 

OICY  J,  T.  ’AA'tt  INGLEY 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENT^v 
SOODARD  3LVD  KING  OF  PRUSSIA 
P.  O.  BOX  8555 

PHILADELPHIA,  PA  19131 

OiCY  ATTN  M.  h.  80RrN£R  SPACE  SCI  -AS 

GENERAL  ELECTRIC  COMPANY 
P.  0.  BOX  1122 
SYRACUSE,  NY  13201 

oicy  at^x  f.  reibert 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  St-jDIES 
816  StAT£  STRE£T  CP.O.  DRAWER  QQ) 

SANTA  3ARBARA,  CA  93102 
OICY  Attn  DA5IAC 
oicy  at^n  don  chandler 

OICY  ATrN  TOM  3ARR£T- 
OICY  At-n  TIM  $tEPhANS 
OICY  ArTN  BARREN  5.  <NAPD 
OICY  AT^N  WILLIAM  MCNAMARA 
OICY  Attn  s.  GAMBIlL 
OICY  attn  MACK  STANTON 

GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 

HM6S 

COURT  STREET 
SYRACUSE,  NY  13201 

OICY  attn  g.  millman 

GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.  0.  90X  6770 
SANTA  3AR8ARA,  CA  93111 

OICY  AttN  JOHN  ISE  jR 
OICY  Att\  joEl  GARBARINO 


geophysical  ;NS^ruTE 

UNIVERSE  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  AT^n:  SECURITY  OFFICER) 
OICY  ATtN  N.  DAVIS  (uncl  only) 

OICY  NEAL  3ROWN  (UNCL  ONLY) 

OICY  at-n  technical  LIBRARY 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EAStf*n  DSv 
77  A  STREET 
NEEDHAM,  MA  02194 

OICY  AT**  MARSHAL  CROSS 

ILLINOIS,  UNIVERSITY  OF 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
URBANA,  IL  61805 

OICY  AT^N  *  r£H 

Illinois,  umversity  of 
107  COBLE  hall 
301  S.  WRIGHT  STREET 
UR 9 ANA,  II  60680 

(All  corres  attn  security  supervisor  cor) 

OICY  ATtN  K.  TEH 
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institute  for  defense  analyses 

400  ARMY -NAVY  DRIVE 
ARLINGTON,  VA  22202 

OICY  ATTN  J.  M.  AElN 
OiCT  ATTN  ERnE5t  BAUER 
OICY  ATTN  HANS  'WOLF HARD 
OICY  ATTN  JOEL  8ENGSTON 

HSS,  INC . 

2  ALFRED  CIRCLE 
BEDFORD,  HA  01730 

OICY  AttN  DONALD  HANSEN 


INTL  TEL  &  TELEGRAPH  CORPORATION 
500  wASHINGTON  AVENUE 
NUNLEY,  NJ  07110 

OICY  ATrN  TECHNICAL  LIBRARY 


jaycoR 

1401  CAM I. NO  DEL  HAR 
DEL  MAR,  CA  92014 

OICY  ATTN  S.  R.  GOLDMAN 


JOHNS  HOPKINS  UN l VERS l TY 
APPLIED  PHYSICS  LABORAtQRY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MO  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 

oicy  at^n  thomas  po?cmRA 

OICY  ArrN  oOHN  DASSOULAS 

LOCKHEED  MISSILES  8  SPACE  CO  INC 
p.  O.  BOX  5Q4 
SUNNYVALE,  CA  94088 

OICY  A^n  DEPT  60-12 
OICY  AT^N  D.  R.  CHURCH  ILL 

LOCKHEED  MISSILES  AND  SPACE  CO  INC 
3251  STREET 

PALO  ALTO,  CA  94304 

OICY  AT^N  MARTIN  WALT  DEPT  52-10 
OICY  ATTn  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTN  W «  L*  iMHOF  DEPT  52-12 

<AMAN  SCIENCES  CORP 
P.  O.  BOX  7463 
COLORAOO  SPRI.NGS,  CO  80953 
OICY  ATtn  t.  MEAGHER 

LlMCABlT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

OICY  ArvN  IRWIN  JACOBS 

M.I.T.  LINCOLN  LABORATORY 
P.  O.  BOX  ?3 
^EXlN&TQN,  ,HA  02173 

OICY  AT^n  DAVID  M.  TOWLE 
OICY  A^^n  P.  hALDRON 
OICY  AT^n  L.  LOUGHLIN 
OICY  A'rTN  0.  CLARK 

MART MARIETTA  CORP 
ORLANDO  DIVISION 
P.  0.  BOX  5837 
ORLANOO,  pL  52805 

OICY  A^N  R.  HEFFNER 

MCDONNEuL  DOUGLAS  CORPORA T I ON 
5501  SOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULE 
OICT  Attn  GEORGE  MROZ 
OICY  AT^n  W.  OLSON 
OICY  At?n  R.  W.  halprin 
OiCY  AT^N  TECHNICAL  library  services 


MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

Santa  Barbara,  ca  93101 


OICY 

ATTN  P.  FISCHER 

Q1CY 

attn  w.  p.  crevier 

OICY 

ArTN  STEVEN  L.  GlTSCHE 

OICY 

0.  5APPENF  lELO 

OICY 

ATTN  R.  BOGuSC-H 

OICY 

A^n  R.  hENDRICK 

OICY 

Attn  RALPH  <l.B 

OICY 

AT^  DAVE  SOWLE 

OICY 

ATTN  F .  FA JEN 

OICY 

AttN  h.  SCHEI3E 

OICY 

A^N  CONRAD  L.  JONG M[RE 

OICY 

ATtn  BARREN  A.  SCHLU£t£r 

MI rRE  CORPORATION,  THE 

P.  0,  30X 
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BEDFORD,  MA  01730 

OICY 

Attn  JOHN  MORGANS TERN 

OICY 

AT^n  G.  HARDING 

OICY 

attn  C.  E.  Callahan 

mitre  CORP 

wESTGATE  RESEARCH  PARK 

1820  DOllY 

MADISON  BLVD 

MCLEAN,  VA 

22101 

OICY 

Attn  w.  hall 

OICY 

Attn  w.  cOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERF I ELD  BLVD. 

SANTA  MONICA,  CA  90404 

OICY  Attn  E.  C.  FIElD  ^R 

PENNSYLVANIA  STATE  JNlVt:RSITY 
IONOSPHERE  RESEARCH  lAB 
313  ELECTRICAL  ENGINEERING  EAST 
UNlVERStTY  PARK,  PA  16802 

(NO  CLASSIFIED  *0  ^IS  ADDRESS) 

OICY  AT^n  IONOSPHERIC  RESEARCH  l.AB 

photqmetrics,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

OICY  Attn  IRVING  L.  KOFSKT 

PHYS l CAL  DYNAMICS  INC. 

P.  0.  BOX  3027 
BELLEVUE,  wA  98009 

OICY  at-hm  E.  0.  FREMOUW 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  10367 
OAKLAND,  CA,  9^610 

ATtn:  A,  rHCMSON 


HD  ASSOCIATES 


p.  0.  BOX 
MARINA  DEL 

9695 

REY, 

ca  90291 

OICY 

AT^ 

F0RRE$t  GILMORE 

OICY 

attn 

BRYAN  GABBARD 

OICY 

A^ 

WILLIAM  5.  wRIGHT  vjR 

Oicr 

attn 

ROBERT  p  _ElEVIER 

OICY 

attn 

WILLIAM  J.  KARZAS 

01CV 

attn 

H.  DRY 

OICY 

attn 

C.  MACDONALD 

OICY 

attn 

R .  TURCO 

RANO  CORPORATION,  the 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

OICY  Attn  CULLEN  CRAIN 
OICY  Attn  ED  9EDROZ I  AN 

RIVERSIDE  RESEARCH  INS^ITc 
80  «EST  END  AVENUE 
NEW  YORK,  NY  10025 

OICY  At^n  VINCE  ’"RAP AN  I 
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SCIENCE  APPLICATIONS,  INC. 

P.  0.  SOX  2351 
LA  JOLLA,  CA  92038 

01CY  Ar-N  LEV*  IS  M.  L  l  NS  ON 
OICY  ATTN  OANIEL  A.  HAMLIN 
01CY  ATTN  0.  SACHS 
01CY  ATTN  E.  A.  STRAXER 
01CY  ATTN  CURTIS  A.  SMITH 
01 CY  ATTN  JACK  MCDOUGALL 

RAYTHEON  CO. 

528  80S T ON  POS'r  ROAO 
SU06URY,  MA  01776 

01CY  AtTn  8AR8ARA  AOAMS 

Science  Applications,  Incorporated 
1710  Goodridge  Drive 
McLean,  VA  22102 

Attn:  J.  Cockayne 

Lockheed  Missile  &  Space  Co., Inc. 
Huntsville  Research  &  Engr .  Ctr. 
4800  Bradford  Drive 
Huntsville,  Alabama  35807 
Attn:  Dale  H.  Davis 


SRI  INTERNATIONAL 
535  RAvEnSwOOO  avenue 


MENLO  PARK, 

,  CA 

94025 

OICY 

ATTN 

DONALD  NEtLSON 

01CT 

AT^n 

AL AN  BURNS 

OICY 

Attn 

G.  SMITH 

OICY 

A^tn 

L.  L.  COBB 

OICY 

at^n 

DAVID  A.  JOHNSON 

OICY 

attn 

WALTER  G.  ChESNUT 

OICY 

attn 

CHARLES  L.  RINO 

OICY 

attn 

WALtER  JAYE 

OICY 

at^n 

M.  BARON 

oicr 

AttN 

RAY  L.  LEADABRAND 

OICY 

A^n 

G.  CARPENTER 

OICY 

ATTn 

G.  PRICE 

OICY 

attn 

J.  PETERSON 

OICY 

ATTN 

R.  HAKE,  JR. 

OICY 

attn 

V.  GONZALES 

OICY 

attn 

D.  MCDANIEL 

TECHNOLOGY 

INTERNATIONAL  corp 

75  WIGGINS 

avenue 

BEDFORD,  MA  01730 

OICY 

attn 

W.  P.  BOQUIST 

TRW  DEFENSE  &  SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONOO  BEACH,  CA  90278 

01CY  ATrN  R.  K .  PLE6UCH 
01CY  AT^n  S.  4ltSCHULER 
01CY  A^N  d.  dee 


VI S I DYNE,  INC. 

19  third  avenue 

NORrH  WEST  INDUSTRIAL  PARK 
BURLINGTON,  MA  01803 

oicy  attn  Charles  Humphrey 
oicy  attn  j.  w.  carpenter 
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IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
UNCLASSIFIED  ONLY 


PLEASE  DISTRIBUTE  ONE  COPY  T0  EACH  OF  The  FOLLOWING  PEOPLE: 

ADVANCED  RESEARCH  PROJECTS  AGENCY  (ARPA) 

STRATEGIC  TECHNOLOGY  OFFICE 
ARLINGTON,  VIRGINIA 

CAPT.  DONALD  M.  LEVINE 


COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  the  NAVY 
WASHINGTON,  D.C .  20560 

DR.  T.  C2UBA 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

DR.  P.  MANGE 
DR.  R.  MEIER 

DR.  £ .  SZUSZCZEWIC2  -  COOE  4127 

DR.  TIMOTHY  COFFEY  -  CODE  4700  (25  COPIES) 

DR.  S.  OSSAJCOW  -  COOE  4780  (100  COPIES) 

DR.  J.  GOODMAN  -  CODE  7560 


SCIENCE  APPLICATIONS,  INC. 
1250  PROSPECT  PLAZA 
LA  JOLLA,  CALIFORNIA  92037 

DR.  D.  A.  HAMLIN 
DR.  L.  LINSON 
DR.  D.  SACHS 


DIRECTOR  OF  SPACE  AMD  ENVIRONMENTAL  LABORATORY 
NOAA 

BOULDER,  COLORADO  80302 

DR.  A.  GLENN  JEAN 
DR.  G.  W.  ADAMS 
DR.  D.  N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 


HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MASS.  02138 

DR.  M.  B.  MCELROY 
DR.  R.  L INOZEN 


PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PENNSYLVANIA  16802 

DR.  J.  S.  NIS6ET 
DR.  P.  R.  RCHR BAUGH 
DR.  D.  E.  BARAN 
DR.  L.  A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BEMSIETT 
DR.  E.  KLEVANS 


UNIVERSITY  OF  CALIFORNIA,  LOS  ANGElES 

405  HULGARD  AVENUE 

LOS  ANGELES,  CALIFORNIA  90024 

DR.  F.  V.  CORONIT l 
DR.  C.  KEWEl 


A.  F.  GEOPHYSICS  LABORATORY 
L.  G.  HANSOM  FIELD 
BEDFORD,  MASS.  01730 


UNIVERSITY  OF  CALIFORNIA,  BERKELEY 
BERKELEY,  CALIFORNIA  94720 

DR.  M.  HUDSON 


DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.  M.  FORBES 
DR.  T.  J.  KENESHEA 
DR.  J.  AARONS 


OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VIRGINIA  22217 

DR.  H.  MJLLANEY 


UTAH  STATE  UNIVERSITY 
4TH  N.  AND  8Th  STREEtS 
LOGAN,  UTAH  84322 

DR.  P.  M.  BAMCS 
DR.  R.  HARRIS 
DR.  V.  PETERSON 
DR.  R.  MEGILL 
DR.  K.  BAKER 


CORNELL  UNIVERSITY 
ITHACA,  NEW  YORK  14850 


COMMANDER 

NAVAL  ELECTRONICS  LABORATORY  CENTER 
SAN  DIEGO,  CALIFORNIA  92152 


DR.  W.  E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 


NASA 

WR.  R.  ROSE  GODDARD  SPACE  FLIGHT  CENTER 

GREENBELT,  MARYLAND  20771 

DR.  S.  CHANDRA 
DR.  K.  MAEDO 


DR.  J.  HEINERL 


U.  S.  ARMY  ABERDEEN  RESEARCH  AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MARYLAND 
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PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  N£w  JERSEY  08540 

OR.  F.  PE  WINS 
OR.  e.  FRIEMAN 


Institute  f-OR  DEFENSE  ANALYSIS 
400  ARMY  NAvY  DRIVE 
ARLINGTON,  VIRGINIA  22202 

DR.  c.  QAuER 


UNIVERSITY  OF  MARYLAND 
COLLEGE  PARn,  MO  20742 
DR.  paPADOPCHJLOS 
DR.  E .  OTT 

UNIVERSITY  of  PITTSBURGH 
PI TTSBGRGn,  =A.  IS2L3 

DR.  N.  CARuS<y 

or.  m.  aiofci 

JEFTv'it  JUC  JME N T  A  1 1 JN  CENTER 
CAMER.A  j  f  A  T I  A 
Alexandria,  va.  :;mu 

i:  cjpil:  if  open  publication 

OTHERWISE  :  COPIES'  I2CY  ATTN  TC 

UNIVERSITY  OF  CALIFORNIA 

lOS  ALAMOS  SCIENTIFIC  LABORATORY 

J-IU,  MS-^o4 

lOS  alamos.  N C*  ME  x ICO  0  7545 
M .  pjngratz 

0.  SIMONS 
s.  sJARASCH 
L  .  DUNCAN 

University  of  California, 

San  Diego 

Dept,  of  Electrical  Engineering 
&  Computer  Sciences 
Mail  Code  C-014 
La  Jolla,  CA  92093 


Dr.  Henry  G.  Booker 
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